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1.  Physical  Properties  of  Fluid  Hydrogen 

1. 1  Heat  Capacity  and  PVT  Measurements 

Dr.  R.  D.  Goodwin  (Project  Leader),  D.  E.  Diller, 

H.  M.  Roder,  Dr.  L.  A.  Weber,  Dr.  B.  A.  Youn glove 

The  PVT  measurements  on  parahydrogen  are  being  conducted 
by  D.  E.  Diller  and  L.  A.  Weber  with  assistance  from  R.D.  Goodwin 
and  occasionally  from  H.  M.  Roder.  H.  M.  Roder  is  programming 
and  conducting  the  computations  with  the  electronic  machines  and  is 
also  examining  the  data  for  self-consistency.  The  heat  capacity 
calorimeter  design  details  are  being  attacked  by  B.  A.  Youn  glove 
full-time  and  L.  A.  Weber  part-time,  and  the  methods  required  to 
obtain  thermodynamic  functions  from  experimental  data  have  been 
surveyed  by  R.  D.  Goodwin. 

Among  the  Appendices  is  a  confirmation,  by  vapor -pres sure 
observations,  that  we  have  indeed  had  parahydrogen  in  the  pipet. 

PVT  data  in  the  range  28  through  70  cc/g.  mol  is  presented  (experi¬ 
mental  runs  Nos.  3  through  31).  Experiments  to  much  higher  densi¬ 
ties  have  been  completed  (through  run  No.  42).  The  highest  densities, 
above  the  point  of  solidification  near  200  atm. ,  are  about  to  be 
measured.  Following  that,  the  low-pressure  measurements  required 
for  computation  of  thermodynamic  functions  will  be  undertaken  in  the 
range  of  amagat  density  from  300  to  30. 

This  portion  of  this  report  was  originally  prepared  as  part 
of  a  different  report  series.  References  to  earlier  reports  in  this 
other  series  occur  in  several  places  in  the  Appendices  and  should 
be  ignored. 
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APPEND1X  I 

Thermodynamic  Functions  from  PVT  Data 

1.  Methods  are  discussed  in  general  and  in  detail  by  authors  of 
references  [1-7]  in  that  order.  The  change  in  thermodynamic  func¬ 
tions  with  pressure  or  density  at  constant  temperature  may  be  com¬ 
puted  from  PVT  data.  The  paths  of  such  computations  are  vertical 
straight  lines  on  figures  1  and  2,  prepared  for  this  discussion.  It  is  * 

necessary  to  choose  some  standard  or  reference  state.  This  leads 
to  confusion  because  in  the  natural  reference  state  of  zero  pressure 
the  functions  A,  S,  and  G  have  logarithmic  infinites.  The  figures 
show  that  isothermal  access  to  the  liquid  region  from  low  densities 
is  blocked  by  the  coexistence  region.  The  choice  of  density,  rather 
than  pressure,  as  an  independent  variable  together  with  temperature 
yields  more  regularly-behaved  functions  for  computational  purposes. 

A  separate  task,  therefore,  is  a  final  interpolation  from  density  to 
pressure  as  independent  variable  for  engineering  applications.  The 
selection  of  functions  and  procedures  for  computation  is  discussed 
below. 

Z.  Reference  or  Standard  States.  Professor  Michels  takes  a  value 
of  zero  for  energy  and  entropy,  E  and  S,  under  standard  conditions, 

0°C  at  1  atm.  To  obtain  values  for  S  along  the  reference  line  of 
amagat  density  p  =  1  at  different  temperatures,  the  specific  heat  as 
a  function  of  temperature  at  p  =  1  is  required  [2].  H.  W.  Woolley 
on  the  other  hand  tabulates  thermodynamic  functions  "for  a  fictitious 
ideal  gas"  at  one  atmosphere  pressure  as  a  function  of  temperature 
[4]  fr om  computations  of  spectroscopic  data  by  statistical  mechanical 
methods  [7,8].  The  situation  of  these  reference  or  standard  states, 
p  =  1  and  P  =  1  atm.  ,  relative  to  the  phases  of  hydrogen  is  shown  by 
the  figures  mentioned. 


-5- 


3.  Basic  Functions.  Residual  functions  are  a  difference  between 
values  for  real  and  ideal  gas  at  the  same  temperature  and  density. 
Whether  explicitly  stated  [  3]  or  not,  they  are  used  in  actual  compu¬ 
tation  [4].  A  defined  residual  function  of  the  data,  Q*,  simplifies 
computations  utilizing  density  as  an  independent  variable  [  3], 

Q*  s  (PV  -  RT)v  =  RT(Z  -  l)v.  (1) 

This  function  is  related  to  isochores,  showing  the  same  regular 
behavior  with  temperature  and  density.  Extrapolation  of  Q*  iso¬ 
therms  to  zero  density  gives  RT*B  where  B  is  the  second  virial 
coefficient.  All  of  the  residual  thermodynamic  functions,  equation 
(13-a),  approach  zero  with  density  (no  logarithmic  infinities).  From 
only  two  such  functions  the  others  may  be  obtained  rather  directly. 
With  density  as  a  variable,  the  residual  Helmholtz  free  energy  may 
be  obtained  without  differentiation  of  the  data, 

u 

A*(u,T)  =  j  Q*dux,  (2) 

o 

where  u  5  1/  v  is  density  and  data  must  be  used  to  sufficiently  low 
density  for  acceptable  error  in  this  residual. 

Any  additional  thermodynamic  function  (variables  u  and  T)now 
requires  differentiation  of  the  data.  It  is  pointed  out  by  F.  Din  [1]  , 
however,  that  it  may  be  preferable  with  precise  and  closely-spaced 
data  to  perform  differentiation  after  an  operation  such  as  equation  (2). 
The  alternatives  are 

%\x 

S*(u,T)  =  (8Q*/8T)udux  (3-a) 

o 

and 

S(u,  T)  =  -  (8A/8T)u. 


(3-b) 
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The  remaining  residual  functions  are 

E*  s  A*  +  TS*  (4) 

H*  =  E*  +  Pv  -  RT  (5) 

G*  =  H*  -  TS*  (6) 

In  (f/P)  =  A*/  RT  -  InZ  +  Z  -  1  (7) 

A  check  is  given  by  application  of 

{8G/8P)t=v  (8) 

to  results  with  pressure  as  independent  variable. 

4.  Specific  Heats.  These  are  highly  exacting  of  data  accuracy  and 


precision: 


C  =  -  T  l  (82Q*/8T2)  du  , 
v  J  u  T 


(9-a) 


C*  -  C*  =  -R  +  T(8P/9T )*/  (8P/8v)t. 


(10-a) 


As  an  alternative,  however,  the  second  differentiation  may  occur  at  a 
late  stage  of  the  computations , 

C  =  (8E/8T)  ,  (9-b) 


C  =  (8E/8T)  , 
v  v 

C  =  (8H/9T)  . 
P  P 


(10-b) 


5.  Residuals  from  Equation  of  State.  When  graphical  methods  are 
employed,  it  is  imperative  for  accuracy  to  handle  only  suitable 
residual  quantities.  At  very  high  densities  any  real  fluid  shows  such 
extreme  deviation  from  ideal  gas  behavior  that  the  conventional 
residuals  [  3]  no  longer  are  relatively  small.  It  is  necessary  then  to 
utilize  an  equation  of  state  for  reference  rather  than  ideal  gas  [6]. 
We  eschew  this  method  since,  for  the  extreme  range  of  densities  to 
be  handled,  it  may  be  expected  that  the  residuals  from  any  possible 
equation  of  state  would  not  be  well-behaved.  Computational  labor, 
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moreover,  becomes  excessive.  The  close-spacing  of  our  data, 
together  with  electronic  computer-methods,  may  permit  use  of  the 
conventional  residuals. 


b.  Michels*  Residuals.  Take  entropy  and  energy  as  zero  at  the 
standard  state,  0®C  and  Pq  =  1  atm.  ,  amagat  density  p  £  v°/v  = 
Entropy  at  any  other  temperature  on  the  p  =  1  isochore  is 


S(l,  T)  ~  S(l, T)  '  S(l, 273) 


•! 


(Cv  dlnT. 

p  =i 

73 


(ID 


Corresponding  Helmholtz  energy  could  be  obtained  with 

(8A/8T)v  =  -S.  (12) 

Professor  Michels  prefers  to  determine  E  on  this  isochore  via  a 
path  through  p  =  0,  [  2].  One  then  has  A  £  E  -  TS. 

The  residual  A*  of  any  thermodynamic  function  A  is  defined 
by 

u 

A*(u,T)  £  ^  [  (8A/8u)t  -  (8A/8u)°  ]  duT,  (13-a) 
o 

where  the  superscript  o  indicates  perfect  gas.  If  relative  values  for 
A  have  been  determined  at  all  temperatures  along  the  p  =  1  isochore, 

values  at  other  densities  on  an  isotherm  are  given  by 

P 

AA  £  A(p ,  T)  -  A(l, T)  =  ^  (8A/8p)TdpT  (14-a) 

i 

P  P 

AA  =  ^  (8A/8p)°  dpT  +  ^  [  (8A/8p)T  -  (8A/8p)°  ]dpT  (14-b) 

l  i 

P 

AA  =  ^  (8A/8p)^dpT  +  A*(p, T)  -  A*(1,T). 

1 


(14-c) 


8- 


If  A  be  Helmholtz  energy,  then 


(8A/8p)°  =  -RT/p 

(15) 

and 

P 

A*(p,T)  =u°^Q*dpT, 
o 

(16) 

where  u°  a  l/v°  is  density  at  standard  state,  0°C,  P  =1  atm. 
o  o  o 

where  molal  volume  is  v°.  Equation(14-c)  then  becomes 

o 

9 

AA  =  -RT  lnp  +  u°  jJa*dpT  “  ^  Q*dpTJ . 

(17«a) 

rp 

A(p,  T)  =  A(l.T)  -  RTlnp  +  u°  \  Q*dpx> 

1 

(17-b) 

For  entropy,  the  relations  corresponding  to  (15)  and  (16)  are 

(8S/8p)°  *  -R/p 

(18) 

and 

P 

S*(p.T)  =  -u°  ^(8Q*/8T)pdpx 
o 

(19) 

wherewith  equation  (14-c)  becomes 

P 

S(p»T)  =  S(1,T)  -  R  In p  -  u°^(8Q*/8T)pdpT.  (20) 

i 

7.  Woolley1  8  Procedure.  This  differs  from  Michels'  procedure  in 

two  essentials:  the  reference  states  are  on  an  isobar  P  =1  atm.  ;  v 

o 

these  states  are  for  hypothetical  ideal  gas,  rather  than  real  gas.  The 
residual  definition  of  equation  (13-a)  may  be  rewritten 

A(u,T)  h  A°(u, T)  +  A*(u,T), 


(13-b) 
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where  A°  is  any  thermodynamic  function  for  perfect  gas.  To  avoid 
the  logarithmic  infinities  in  the  functions  A°,  S°  and  G°  for  the 
perfect  gas  at  zero  pressure,  the  spectroscopically  computed  values 
A°(uo,  T)  at  Pq  =  1  atm.  ,  corresponding  to  uq(T),  are  employed  for 
reference  at  each  T, 

u 

A°(u,T)  =  A°<uo,T)  +  J  (8A/du)°duT.  (21) 

uo 

Equation  (13-b)  then  becomes 

A(u, T)  =  A°(uo.T)  +  J  (8A/8u)°duT  +  A*(u,T),  (22-a) 

uo 

or 

A(p,T)  =  A°(po,T)  +  ^  (8A/8p)°dpT+  A*(p,T),  (22-b) 

Po 

where  pQ,  the  amagat  density  of  perfect  gas  at  Pq  =  1  is 
temperature -dependent. 

If  A  be  Helmholtz  energy,  then  by  equation  (15) 

P 

y  (8A/8p)°dpT  =  -RT  lnp  +  RT  lnpQ.  (23) 

Po 

But  we  find  the  temperature -dependence  of  pQ  as  follows  for  perfect 

gas.  ( 

p  =  v  /v  where  v  =  RT /  P  ,  (24) 

ooo  o  o 

whence 

Po  =  >Povo/RT0»To/T)-  <25> 

Combining  equations  (16),  (23),  and  (25)  with  (22-b)  for  Helmholtz 
energy. 
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A(p,T)  =  A°(po,T)  -  RTlnp  +  RT  [  In  (Pqv°/RTo)  +  ln(TQ/T)] 

P 

+  ^  Q*dpT.  (26) 

o 

The  corresponding  expression  for  entropy  is  derived  by  exact  analogy, 

S<p,T)  =S°(po,T)  -  R In p  +  R[ln(pov°/RTQ)  +  ln(TQ/T)] 

P 

-  u°  f  (8Q*/8T)  dp_.  (27) 

o  J  P  T 

o 

The  first  terms  on  the  right  of  equations  (26)  and  (27)  represent  the 
tabulated  reference  values  for  hypothetical  ideal  gas  at  one  atmos¬ 
phere  pressure. 

8.  Orthobaric  Densities.  Some  curvature  of  isochores  near  satura¬ 
tion,  not  shown  by  figure  1,  may  be  suspected.  If  experimental  iso¬ 
chores  were  measured  with  very  small  temperature  increments 
extending  into  the  coexistence  region,  a  set  of  isochores  would  define 
the  vapor  pressure  curve.  Such  desirable  measurements  usually  are 
not  undertaken  as  a  part  of  PVT  work  because  of  the  greatly  increased 
labor.  If  heats  of  vaporization  and  the  vapor  pressure  curve  (curva¬ 
ture)  are  known,  the  difference  between  orthobaric  densities  may  be 
computed.  Hence  the  density  of  saturated  liquid,  from  triple  to 
critical  points,  is  valuable  data.  This  could  be  utilized  more  directly 
with  the  rule  of  rectilinear  diameter,  provided  two  vapor  densities 
were  known  [15].  Some  such  information  is  required  to  extrapolate 
ordinary  PVT  isochores  onto  the  coexistence  boundary. 

9.  Compressed  Liquid  Region.  Entry  to  this  region  requires  compu¬ 
tation  of  functions  of  density  along  an  isotherm  which  crosses  the 
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coexistence  region.  Applicable  relations  for  the  crossing  are 


*S  S  Q/T, 

(28) 

(8S/8v)t  *  (8P/8T)v, 

(29) 

(8A/8v)t  =  -P. 

(30) 

Upon  integration  of  (29)  and  (30)  in  the  coexistence  region  and  elimi¬ 
nation  of 

1AT  = -Q/(81nP/81nT).  (31) 

Thermodynamic  properties  in  the  compressed  liquid  region  cannot  be 
related  to  any  low-pressure  standard  state  by  means  of  PVT  data 
alone.  Latent  heats  [10,11]  and  the  vapor  pressure  curve  [14]  must 
be  known.  An  alternative  is  to  employ  a  specific  heat  measurement, 
Cv>  extending  from  the  triple  point  through  the  compressed  liquid 
region  to  a  temperature  well  above  critical  where  the  functions  are 
determined  by  PVT  data  alone,  figure  1,[  7]. 

10.  Low  Pressure  Data  Requirements.  While  our  interest  and 
measurements  are  concerned  with  high  density  regions,  we  find  that 
a  self-consistent  table  of  thermodynamic  properties  preferably  is 
based  upon  some  low-density  standard  or  reference  states:  equations 
(26)  and  (27),  for  example.  Just  how  low  in  density  must  experi¬ 
mental  measurements  be  extended?  Isothermal  plots  of  Q*  vs.  p 
become  linear  at  low  density.  The  linear  portions  provide  values  for 
second  and  third  virial  coefficients.  In  principle  these  isotherms 
should  cover  the  temperature  range  of  the  subsequent  calculations. 

In  practice  some  extrapolation  of  the  virial  coefficients  to  lower 
temperatures  may  be  necessary.  Figure  8  of  the  compendium  of 
Woolley  etal.  [4]  shows  that  the  hydrogen  plots  are  linear  in  the 
range  p  <  200.  Our  figure  2  shows  that  this  covers  most  of  the  gas 
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and  vapor  region  of  densities  less  than  critical  (pc«  335).  Our 
apparatus  is  capable  of  measurements  to  densities  as  low  as  p  =  20. 
We  therefore  should  be  able  to  determine  second  and  ttiird  virial 
coefficients  by  means  of  data  from  p  =  20  to  p  =  200,  with  resulting 
isothermal  equations  of  state  valid  for  0  <  p  <  200  excluding  the 
condensed  phases. 

11.  Discussion.  A  commendable  procedure  of  Woolley  etal.  [4]  is  to 
maintain  separate  tabulations  of  terms  in  equations  such  as  (26).  Any 
errors  in  individual  terms  will  not  invalidate  tables  of  the  others. 

Our  figures  show  that  the  p  =  1  isochore  does  not  penetrate  the 
coexistence  region.  For  simplicity  of  conception,  therefore,  the 
virial  equation  could  be  employed  with  equations  (26)  and  (27), 
together  with  spectroscopic  data  for  ideal  parahydrogen,  to  compute 

a  table  of  thermodynamic  function  values  at  integral  tem”. :  -atures 
along  the  p  :1  isochore.  Equations  (17- b)  and  (20)  then  would  be 
applicable  with  the  virial  equation  for  isothermal  computations  as  a 
function  of  increasing  density,  with  upper  limit  the  coexistence 
boundary  or  limit  of  validity  of  the  virial  equation  at  about  p  =  200. 

12.  Tentative  Computational  Tables.  The  following  list  of  tables  may 
serve  as  a  guide  to  the  computational  tasks  required: 

1.  Functions  for  ideal  gas  at  1  atm.  and  integral  tempera¬ 
tures. 

2.  Functions  for  real  gas  at  p  =  1  and  integral  tempera¬ 
tures. 

3.  The  change  in  fucntions  between  p  =  1  and  p  =  p^, 
where  the  integral  values  of  p^  extend  to  the 
coexistence  boundary  for  integral  temperatures 
up  to  critical,  and  to  about  p  =  200  at  higher 
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integral  temperatures,  employing  the  virial  equation 

up  to  "third"  coefficient,  and  direct  PVT  data  for  the 

vapor  region  200  <  p  <  p  . 

c 

4.  The  change  in  functions  upon  crossing  the  coexistence 
region,  vapor  to  liquid,  at  integral  temperatures  up 
to  critical. 

5.  The  change  in  functions  of  compressed  liquid  from 
each  state  of  saturation  up  to  integral  densities,  p^, 
extending  to  the  solid  or  limit  of  PVT  data,  at 
integral  temperatures  up  to  critical. 

6.  The  change  in  functions  of  gas  and  fluid  from  about 
p  =  200  to  integral  values  p^  extending  to  the  limit 
of  PVT  data  at  integral  temperatures  above  critical. 

R.  D.G. 
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Figure  I* 
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Figur*  2- 
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APPENDIX  H 

Expansion  of  Approximate  Equation  of  State 

In  seeking  regularly-behaved  isotherms  for  examination  of  the 
precision  of  experimental  densities,  and  therefore  also  comparison 
with  the  virial  equation,  it  might  be  useful  to  expand  the  equation  of 
state  (Report  13,  page  18)  units  of  cc.  ,  g.mol.  ,  atm.  ,  degree  Kelvin, 

Pv/RT  =  v/(v  -  b  +  c/v)  -  (A/RTvm-1)  [1  -  (vQ/v)n].  (1) 

b  =  26.  9  vQ  =  23.478 

c  =  237.  5  m  =  1.86 

A  =  125,000  n  =6 

The  expansion  of  the  first  term  of  (1)  converges  for  molal  volumes 
greater  than  7  cc.  , 

3  4  5 

v/(v  -  b  +  c/v)  =  1  +  b,u  +  cxu2  +  d,u  +  eju  +  f,u  +  ...,  (2) 

where  the  reduced  density  is 

u  s  v  /v;  v  =  66.  95  cc.  , 
c  c 

and  the  constants  are  bt  =0.  4018, 

c,  =  0.1085  e,  =  0.  00321 

dj  =  0.0223  fj  =  0.  00011. 

Reduced  constants  for  the  second  term  of  (1)  are 

k,  sA/Rv™-1  =  40.987  k,k2  =  0.07621. 

k2  s  (v  /v  )6  =  0.  001860 

c  o  c 

The  expanded  equation  then  is  of  the  form 

Pv/RT  =  1  +  bjU  +  CjU2  +  •  •  •  -7jfum  I[l-k2u  ].  (3) 

For  approximations  it  is  recognized  that  m  differs  but  slightly  from 
2.  (The  constant  kj  is  a  function  of  m. )  Combining  similar  terms. 
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Pv/RT  =  1  +  (bj  -  k  um“2/T)  u  +  Cju2  +  . . .  +k1k2um+5/T.  (4) 

The  more  useful  form  of  state  and  virial  equations »  respectively! 
then*  are 

(z  -  l)/u  =  (Vj  -  kjUm’2/T)  +  CjU  +  djU2  +  . . . 

+  kjk^^/T.  (5) 

(z  -  l)/u  =  Bj  +C^u  +  DjU2  + .  (6) 

where  the  conventional  virial  coefficients  B^,  C ^ ,  Da .  .  .  are  in 
reduced  form.  No  simple  comparison  between  coefficients  of 
equations  (5)  and  (6)  is  possible  since  the  coefficients  of  the  latter 
must  be  temperature -dependent,  whereas  for  the  former  they  are  not. 

R.D.G. 
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APPENDIX  IU 

Improved  Temperature-Control  Technique 

Improved  temperature  control  of  the  PVT  pipet  has  been 
achieved  in  practice  by  means  of  appreciable  simplification  of  techni¬ 
que  and  required  apparatus.  The  control  thermocouple  between 
refrigerant  tank  and  pipet  is  not  used.  The  bucking  circuit  for  this 
thermocouple  is  not  used.  Instead,  we  are  using  the  platinum 
resistance  thermometer  for  control  as  well  as  measurement.  Output 
signal  of  the  electronic  galvanometer,  used  with  the  microvolt  poten¬ 
tiometer  for  reading  the  platinum  thermometer,  is  connected  directly 
to  the  power  regulator  for  pipet  heating  temperature  control.  When 
it  is  required  to  standardize  the  potentiometer,  the  steady-state 
heating  current  first  is  carefully  adjusted  on  the  power -regulator. 

The  galvanometer  signal  then  is  disconnected  from  the  power  regu¬ 
lator  during  the  potentiometer  standardizing  procedures,  including 
thermometer  current  standardization. 


O.  E.  D. 
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APPENDIX  IV 

Potentiometer  Battery  Stabilization* 

Drift  rate  under  a  few  parts  per  million  per  hour  is  desired  in 
the  main  working  battery  for  a  microvolt  potentiometer.  The  history 
of  our  observations  follows. 

1)  Nickel-cadmium  alkaline  batteries  were  said  to  be  absolutely 
trouble-free.  Our  4-amp.  -hour  pulton  cells  exude  caustic.  Their 
use  at  full  charge,  maintained  by  trickle -charging  (Report  2,  page  2, 
and  Report  11,  figure  1),  is  unsatisfactory. 

2)  Mercury  batteries,  e.g.  ,  Mallory  RM-42R,  while  excellent 
for  1  milliampere  circuits,  are  unsatisfactory  at  the  5  ma.  -level  for 
the  main  potentiometer  battery. 

3)  Nickel-cadmium  cells  used  in  single-  or  in  two-stage  circuits 
as  recommended  by  E.  E.  Watson  (Report  2,  page  2),  merely  as 
stabilizers  across  another  power  source  at  50  percent  of  capacity 
charge  are  quite  unsatisfactory.  Investigation  shows  that  effective 
internal  resistance  of  0.1  ohm  for  a  6-volt  battery  at  one  ampere 
rises  into  the  range  1  to  10  ohms  at  charge -discharge  rates  in  the 
sub-milliampere  level,  accompanied  by  intolerable  hysteresis  effects. 

4)  A  conventional,  100 -amp. -hour ,  low-discharge  type  of 
lead-acid  battery  in  our  figure  3  stabilization  circuit  is  giving  satis¬ 
faction.  The  battery  is  thermally  and  electrostatically  shielded.  The 
two -microfarad  condensers  greatly  reduce  noise  which  otherwise 
appears  in  the  electronic  galvanometer  amplifier  used  with  the  poten¬ 
tiometer  when  the  rectifier  power  source  is  connected  to  the  battery. 

♦See  R.  J.  Corruccini,  Stabilization  of  a  d.  c.  current  source  by 
means  of  a  parallel-connected  battery,  NBS/CEL.  Dab  Note  60-12, 
October  13,  I960. 
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The  ten-turn  helical  rheostat  permits  any  effective  battery  drift  to 
be  adjusted  to  zero.  Our  battery  is  stable  under  a  charging  current 
of  about  0.  7  milliampere.  The  manufacturer  claims  a  self-discharge 
rate  of  less  than  15  percent  per  year  or  1.  7  milliampere. 

R.  D.  G. 
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STABILIZED  BATTERY  FOR 
MICROVOLT  POTENTIOMETER 


FIGURE  3 


POTENTIOMETER 
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APPENDIX  V 

Pressure  Correction  for  Fluid  in  Pipet  Capillary 

From  page  26,  Report  12,  the  number  of  gram-moles  of  hydrogen 
in  our  transition  capillary  is 

N  =  PV  S/RT  . 

X  X  o 

The  dimens  ion,  less  integral,  0  <  S  <  1,  page  43,  Report  12,  will  be 
taken  here  as  unity.  Pressure  due  to  the  capillary  column  then  is 

P  «  l(f3  •  N  M/A, 
x  x 

where  hydrogen  mol.  wt.  is  M  =  2  and  capillary  area  is  A  =  10  cm2 

.3 

The  conversion  factor  shown  is  10  atm. / (g /  cm2).  Relative  error  is 

P  /p  «  10_3V  M/ART  , 
jy  x  o 

reducing,  with  above  values,  to 

P  /P  «  2  V  /RT  =  (16 /T  )  •  10~4, 

X  X  o  o 

where  Tq  is  pipet  temperature.  This  correction  is  not  applied  by  our 
routine  computation,  page  19,  Report  15. 


R.  D.  G. 
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APPENDIX  VI 

Isochore  for  n-Hydrogen,  Run  No.  2 


It  is  convenient  to  compare  this  experimental  result  with  the 
data  of  Stewart  and  Johnson  and  with  our  equation  of  state ,  provided 
the  data  be  interpolated  to  a  true  isochore.  For  data  of  Exhibit  VIII 
in  Report  15,  an  arbitrary  mean  molal  volume  of  30.6  cc.  is  selected. 
The  required  interpolation  of  pressures  to  this  constant  density  is 
obtained  by  linear  interpolation  of  Stewart  and  Johnson' s  data  between 
v  =  30.0  and  v  =  31.0, 

-AP/Av  =  0.475T-1.4.  (a) 

Data  of  S.  and  J.  at  v  =  30.  6  was  computed  by 


PgJ  =  7.19  T  -  171.2.  (b) 

Our  equation  of  state  at  this  density  is 


P  ,  =  7. 1594  T  -  171.  55.  <c) 

calc.  ' 

The  second  column  of  table  1  gives  the  experimental  pressures  cor¬ 
rected  to  a  constant  volume  of  30.6,  P^^  by  means  of  equation  (a). 
A  deviation  plot  of 


(P30.6  "  Pcalc.)/P30.6 

vs.  temperature  shows  that  this  isochore  in  the  compressed  liquid 
region  curves  toward  the  critical  isochore  with  increasing  temperature. 


R.D.G 
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Table  l 

Isochores  for  n-Hydrogen,  Run  No.  2 
at  v  a  30.  6  cc.  /mol,  Appendix  VI 


T°K 

p 

30.6 

28 

29.001 

30 

43. 558 

32 

58.027 

34 

72. 492 

36 

86.  999 

40 

115.878 

45 

151. 365 

50 

186. 369 

55 

220. 996 

60 

255. 180 

65 

288.  967 

PSJ 

Pcalc 

30.  1 

28.  91 

44.  5 

43.23 

58.9 

57.55 

73.  3 

71.87 

87.6 

86. 19 

116.4 

114.83 

152.4 

150.62 

188.  3 

186.42 

(224.  3) 

222.22 

(260.  2) 

258.01 

(296.  2) 

293.81 
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Vapor  Pressure  Observations,  Parahydrogen 

The  purified  hydrogen  used  for  the  PVT  measurements  was  con¬ 
verted  to  the  para  form  by  passing  slowly  over  a  ferric  oxide  catalyst 
at  about  19.  6aK.  The  efficiency  of  the  conversion  was  checked  in 
several  cases  when  the  sample  was  cooled  below  the  point  where  the 
isochore  intersects  the  vapor  pressure  curve.  The  resulting  vapor 
pressures  were  compared  with  values  for  20.  5 °K- equilibrium  hydrogen, 
e-H2,  from  the  literature.  The  table  following  illustrates  the  agree¬ 
ment  of  our  data  with  that  of  Hoge  and  Arnold  [1].  The  corresponding 
vapor  pressures  for  normal  hydrogen,  from  White,  Friedman,  and 
Johnston  [  2]  ,  are  included  for  comparison. 


L.  A.  W. 
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Table  2 

Vapor  Pressure  Observations 


Run 

T,#K 

P(this  work) 
(atm.) 

e-H2[l] 

(atm.) 

n-H2[2] 
(atm. ) 

5 

Zb 

3.  9811 

3.  9777 

3.  8249 

6 

Z1 

4.8215 

4.  8226 

4. 6296 

14 

30 

8. 1041 

8.  1087 

7.  7786 

25 

32 

11.051 

11.033 

10.630 

9 
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APPENDIX  YU! 

PVT  Data  for  Parahydrogen 

Tables  3  through  31  present  experimental  results  in  the  range 
28  through  70  cc/g.mol.  The  following  limitations  should  be  care- 
hilly  considered  relative  to  any  utilization  of  these  data: 

1)  A  small,  indeterminate,  but  not  insignificant  error  exists 
in  the  temperatures,  as  described  below.  Some  of  these 
experiments  therefore  will  be  repeated. 

2)  These  data  have  not  been  smoothed. 

3)  The  initial  point  of  a  few  runs  falls  on  the  vapor  pressure 
curve. 

These  tables  have  been  computed  from  laboratory  "input" 
observations  as  described  by  Exhibits  VI  and  VII  of  Report  15.  This 
laboratory  has  changed  from  an  IBM-650  computing  machine  to  a 
CDC-1604  during  the  course  of  the  experiments.  The  latter  was  used 
to  print  out  the  data  tables.  Greater  flexibility  of  the  latter  computer 
is  reflected  by  our  revised  laboratory  data  sheets.  Exhibits  A  and  B. 

Temperature  Error.  Initial  measurements  purportedly  showed 
absence  of  significant  spurious  potentials  from  the  platinum  thermo¬ 
meter  measuring  circuit  when  the  thermometer  current  was  off.  After 
the  PVT  measurements  were  completed,  a  spurious  residual  potential 
was  found  which  drifted  from  0.  36  to  0.  24  microvolt  in  eight  hours 
under  conditions  used  for  PVT  measurements.  This  potential  was  of 
the  same  sign  as  those  produced  for  measurement  by  the  one  milli- 
ampere  thermometer  current.  Temperatures  given  by  tables  3  through 
31  therefore  are  higher  than  the  true  temperature.  The  following  table 
illustrates  the  magnitude  of  the  error,  which  must  be  regarded  as 
indeterminate. 
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Temperature  Error  Corresponding  to  0.  3  Microvolt 
Residual  Signal  on  Platinum  Thermometer  Measured 
with  One  Milliampere  Current 


T  dR/dT  Error,  *K 


15 

.0081 

0.0370 

20 

.  0186 

0.0161 

25 

.  0329 

0.0091 

30 

.  0483 

0. 0062 

35 

.  0632 

0. 0047 

40 

.  0761 

0.  0039 

60 

.  1045 

0. 0029 

100 

.  1102 

0. 0027 

If  the  temperature  error  be  ignored,  a  corresponding  error  in  pressures 
may  be  considered.  Rough  values  for  the  worst  condition  (occurring 
at  a  low  pressure  of  5  atm.  and  an  isochore  slope  of  10  atm.  per  degree) 
yield  at  about  17°K  an  error  in  the  order  of  6  percent  in  pressure. 

Three  examinations  for  self-consistency  or  mistakes  have 
been  performed. 

1)  The  laboratory  P  vs.  T  data  for  each  run  is  difference  prior 
to  any  PVT  computation  as  in  Exhibit  IV,  Report  15.  A  plot  of  these 
first  differences  in  the  form 

AP/AT  vs.  T*  (Tz  +  T,)/2 

affords  a  sensitivity  of  examination  equal  to  or  greater  than  the  pre¬ 
cision  of  measurement  of  these  variables.  A  few  early  runs  were 
repeated  because  of  inadequate  smoothness  in  these  plots. 

2)  For  each  run,  the  PVT  computation  is  carried  out  to  obtain 
the  nearly  constant  density  at  each  temperature.  Each  of  these  densities 
depends  upon  numerous  small  apparatus  corrections,  such  as  the 
temperature  of  the  external  null-detector  diaphragm.  Any  erroneously- 
recorded  laboratory  data  will  produce  irregularities  in  some  order  of 
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the  differences  of  these  densities  within  a  single  ran.  Exhibit  V, 

Report  51,  detected  usually  by  a  plot  of  Av/&T  vs.  T. 

3)  Self-consistency  of  the  density  between  different  runs  at 
different  densities  is  equivalent  to  such  consistency  in  measurements 
of  the  amount  of  gas  released  from  the  pipet  at  the  end  of  each  run. 
After  the  PVT  computations  for  a  set  of  numerous  runs  are  complete, 
it  is  required  to  plot  and  examine  for  smoothness  a  suitable  type  of 
isotherm  containing  one  point  from  every  run.  To  ask  for  an  isotherm 
linear  in  density  is  to  ask  for  an  equation  of  state.  We  have  found 
adequate  sensitivity  for  the  present  purpose,  despite  the  curvature  at 
high  densities,  by  using  the  value  of  the  second  virial  coefficient,  B, 
for  normal  hydrogen  together  with  our  parahydrogen  data  Z  *  Pv/RT 
and  plotting 

[  (Z  -  l)v  -  Bjv  vs.  1/v. 

The  results  in  tables  3  through  31  were  so  examined  on  the  60 *K 
isotherm. 

H.  M.  R.  and  R.  D.  G. 
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Since  the  pressure  in  the  thermal  conductivity  apparatus 
reaches  300  atmospheres,  it  was  necessary  to  open  the  sealed  plati¬ 
num  resistance  thermometers  to  the  surrounding  atmosphere  to 
prevent  them  from  being  crushed  even  though  this  would  also  cause 
them  to  lose  calibration. 

To  recalibrate  the  thermometers,  the  hot  plate  and  the  cold 
plate,  each  containing  a  platinum  resistance  thermometer,  were 
placed  in  physical  contact.  The  hot  plate,  cold  plate,  and  block  were 
then  in  thermal  equilibrium.  The  measuring  block  could  be  held  at 
any  desired  temperature  and  the  two  thermometers  inside  checked 
against  a  sealed  resistance  thermometer  with  a  previous  calibration 
outside  the  pressure  cell. 

Unfortunately,  we  have  not  received  the  calibration  on  the 
external  thermometer;  therefore,  all  the  data  presented  in  tables  I 
and  II  must  be  relative  and  is  given  in  terms  of  resistance  as  a 
function  of  the  resistance  of  the  calibrated  therr  neter  rather  than 
a  function  of  temperature. 

The  calibration  was  done  for  0.  1,  1.  0,  and  10.0  atmospheres 
with  helium  gas  in  the  system.  The  block  was  brought  to  a  specified 
temperature  and  held  at  that  temperature  while  the  readings  were 
taken;  then  the  pressure  was  changed  and  the  temperature  held 
constant.  This  was  rather  easy  to  accomplish  since  the  heat  capacity 
of  the  block  is  many  times  that  of  the  gas  in  the  system.  The  tempera¬ 
ture  can  be  held  to  at  least  ±0.  001°K  or  ±0.  0001  ohms  on  the  resistance 
thermometer.  Many  times  the  block  was  steady  to  ±0.  0001*K. 
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The  resistances  of  the  platinum  thermometers  were  measured 
by  a  potentiometric  method.  A  current  of  2  milliamperes  ±5  micro¬ 
amperes,  supplied  by  three  mercury  cells  in  series  with  2,000  ohms 
and  measured  by  a  1-ohm  standard  resistor,  was  passed  through  the 
thermometers.  The  apparent  voltage  on  the  platinum  resistance 
thermometer  was  measured  both  with  and  without  current  flowing 
through  it,  and  the  difference  between  these  two  measurements  was 
the  true  voltage  on  the  thermometer.  A  slight  correction  to  the  current 
was  also  necessary  since  the  1-ohm  resistor  was  actually  0.  99998  ohms. 
The  ratio  of  the  corrected  voltage  to  the  corrected  current  then  gave 
the  resistance  of  the  platinum  thermometer. 

Preliminary  graphs  of  the  data  indicate  that  there  is  some  con¬ 
duction  of  heat  from  the  thermometers  by  the  gas.  This  may  mean 
there  will  be  a  slight  error  produced  in  the  temperature  measurement 
which  will  be  a  function  of  the  thermal  conductivity  of  the  gas  in  the 
system.  For  this  reason,  and  also  to  make  various  other  calibration 
checks  on  the  system,  we  will  measure  the  thermal  conductivity  of  air 
first  and  then  check  the  results  of  the  system  against  the  results  of 
other  researchers. 

Fortunately,  the  minor  disturbances,  such  as  thermocouple 
breaks  and  vacuum  failures,  which  have  slowed  the  progress  of  the 
apparatus  are  becoming  less  frequent. 


Table  I 

Platinum  Resistance  Thermometer  Calibrations  -  Hydroge 
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Platinum  Resistance  Thermometer  Calibrations  -  Nitrogen  Range 
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2.  Cryogenic  Instrumentation 

J.  Macinko,  P.  Smelser,  R.  C.  Muhlenhaupt, 

C.  E.  Miller,  Dr.  R.  B.  Jacobs 

The  main  effort  of  the  instrumentation  project  during  the  past 
quarter  has  been  devoted  to  the  design  and  construction  of  test 
apparatus  with  a  secondary  effort  applied  to  the  survey  of  flowmeters. 
Manufacturers  have  been  requested  to  submit  instruments  on  a  con¬ 
signment  basis  for  testing.  Upon  completion  of  testing,  the  instru¬ 
ments  will  be  returned  and  results  and  recommendations  will  be  sent 
to  the  manufacturers. 

Pressure  Transducers 

The  test  apparatus  (figure  11,  First  Quarterly  Report),  has 

» 

been  built,  installed  at  the  test  site  and  the  associated  control  and 
recording  instruments  have  been  checked  out.  The  apparatus  has 
been  cold  shocked  and  mock  tests  have  been  run.  A  delay  in  the 
shipping  of  the  pressure  standards  prevented  any  tests  from  being 
performed  to  date.  However  the  test  gages  have  been  installed  and 
preliminary  tests  are  being  run  while  waiting  for  the  dead  weight 
calibration  tester.  The  overall  performance  will  be  evaluated  after 
calibration. 

Each  transducer  will  be  calibrated  at  temperatures  of  300°  K, 
90°K,  76°K  and  20°K.  In  addition,  the  temperature  bath  can  be 
controlled  in  the  ranges  of  20  to  33° K  and  76  to  100°K.  Temperature 
coefficients  in  these  regions  can  thus  be  obtained.  Temperature 
cycling  and  hysteresis  effects  will  be  determined  in  addition  to 
repeatability,  zero  shift  and  linearity. 

Temperature  Sensors 

A  test  apparatus  (figure  1)  for  determining  the  response  time 
of  temperature  sensors  has  been  designed  and  construction  will  begin 
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Figure  1 

TEMPERATURE  SENSOR  RESPONSE  TESTER 
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shortly.  This  unit  has  a  spring  loaded  plunger  which  will  move  a 
temperature  sensor  from  a  controlled  temperature  in  the  vapor  region 
into  a  controlled  liquid  bath.  This  will  give  the  response  time  for 
sensors  subjected  to  a  vapor  to  liquid  change  in  surroundings.  As 
designed,  the  test  chamber  will  accommodate  sensors  up  to  3/4"  in 
diameter  and  up  to  20"  in  length. 

A  second  unit  is  being  designed  to  utilize  the  rapid 
change  in  liquid  temperature  brought  about  by  adiabatic  compression. 
With  the  sensor  submerged,  a  temperature  step  function  can  be 
obtained  in  the  ambient  liquid  and  the  response  time  determined. 
Several  basic  design  problems  must  be  resolved  before  the  overall 
design  takes  its  final  form.  Two  apparatuses  concerned  with  the 
dynamic  characteristics  of  temperature  sensors  are  being  used  in  an 
effort  to  obtain  data  directly  applicable  to  actual  situations,  and  to 
ascertain  the  feasibility  of  describing  character’.*  Me  s  independent 
from  the  heat  transfer  characteristics  of  the  test  fluids. 

Liquid  Level  Gages 

The  instrumentation  program  will  include  the  calibration  and 
evaluation  of  liquid  level  gages.  Present  plans  call  for  a  primary 
standard  capable  of  accurately  following  a  changing  liquid  level. 
Several  ideas  for  a  standard  are  being  investigated  and  tested.  As 
soon  as  a  decision  is  made  on  the  standard,  design  of  the  test 
apparatus  will  begin. 

Tests  have  just  been  completed  on  a  transistor  liquid  level 
indicator  using  carbon  resistors  as  sensing  elements.  A  schematic  of 
the  circuit  is  given  in  figure  2.  Factors  affecting  choice  of  resistor 
size,  type  of  liquid,  stability,  reliability,  and  other  considerations 
are  discussed  in  a  paper  which  is  undergoing  final  revision  and  will 
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LIQUID  LEVEL  INDICATOR  FOR 
LIQUID  HELIUM,  LIQUID  HYDROGEN,  AND  LIQUID  NITROGEN 
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be  published  shortly  in  "Review  of  Scientific  Instruments".  A  copy 
of  the  paper  will  be  included  in  a  future  report. 

Flowmeter  Survey 

Due  to  the  lack  of  flow  test  facilities  at  present,  and  no 
definite  plans  for  any  future  facility,  only  a  limited  amount  of  test 
data  is  contained  in  the  flowmeter  survey.  The  survey  compiled  to 
date  is  based  on  manufacturers'  literature,  results  of  applications 
of  meters  in  industry,  a  limited  number  of  test  reports,  and  several 
other  surveys  of  flowmeters  available  from  technical  publications. 
Inquiries  have  been  sent  to  most  of  the  manufacturers  of  flowmeters 
applicable  to  cryogenic  fluids;  a  few  more  will  be  sent  out  shortly. 
The  inquiries  have  requested  a  physical  description,  operating 
principles  and  characteristics  including  range  and  accuracy,  and 
price  if  commercially  available.  The  first  draft  of  the  survey 
contains  the  information  available  at  present  and  will  be  distributed 
shortly;  continuous  additions  and  revisions  are  being  made  as  more 
information  is  received. 

2.  1  Forced  Vibration  Densitometer  Studies 

A  test  apparatus  similar  to  the  one  shown  in  figure  12,  First 
Quarterly  Report,  has  been  designed  and  built.  Tests  have  been 
performed  on  the  apparatus  without  fluid  to  establish  the  dynamic 
characteristics  of  the  instrument.  Results  of  these  tests  indicated 
several  undesirable  characteristics  to  be  present.  The  most  critical 
of  these  were  frequent  failures  of  the  brass  bellows  due  to  over¬ 
stressing  and  work  hardening,  and  harmonics  which  cause  distortion 
in  the  dynamometer  output  signal. 

The  following  modifications  were  made  in  an  attempt  to 
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a.  The  brass  bellows  were  replaced  by  bellows  of  beryllium- 

copper. 

b.  The  amplitude  of  the  oscillating  passage  was  decreased 
from  0.  22"  to  0.  12". 

c.  Improved  passage  guides  were  designed  and  installed  to 
insure  only  one  degree  of  freedom  for  the  passage. 

A  Statham  dynamometer  was  selected  to  measure  the  force 

transmitted  to  the  passage.  This  gage  has  a  flat  frequency  response 

4 

up  to  3000  cps,  a  stiffness  of  6  x  10  #/in.  and  a  load  capacity  of  ±10 
lbs.  It  is  assumed  that  with  this  stiffness  the  driver  and  passage 
motion  may  be  considered  in  phase.  A  proximity  pickup  was  installed 
to  check  the  validity  of  this  assumption.  Calibration  tests  will  be  run 
as  soon  as  the  instrument  proves  to  have  sui table  endurance  and 
undesirable  harmonics  have  been  eliminated. 

3.  Cryogenic  Design  Principles  and  Materials  Utilization 
D.  B.  Chelton,  L.  E.  Scott,  J.  A.  Brennan  and 
B.  W.  Birmingham 

General  assistance  on  Project  Centaur  has  continued  during 
the  reporting  period.  Emphasis,  however,  has  been  placed  on  the 
design  of  low  temperature  observation  windows  for  various  containers, 
the  preparation  for  testing  of  bearings  for  the  zero-gravity  centrifugal 
vent  device  and  the  adaptation  of  NBS  developed  seals  for  fuel  line 
flange  application. 

A  detailed  discussion  of  some  of  the  technical  problems  that 
have  been  considered  on  Project  Centaur  are  described  below.  In 
addition  to  those  described,  discussions  have  been  held  with  CVA 
personnel  on  other  matters  of  general  cryogenic  importance  such  as 
properties  of  materials,  temperature  stratification,  and  instrumenta¬ 
tion.  Three  visits  to  Convair  Astronautics  (CVA)  were  made  by  NBS 
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personnel.  Several  telephone  conversations  to  assist  CVA  personnel 
have  taken  place.  Also,  CVA  has  visited  NBS  for  assistance  and  to 
observe  certain  tests. 

3.  1  Observation  Windows 

Early  in  this  reporting  period,  a  CVA  window  frame  for  the 
Zero-G  test  container  was  brought  to  Boulder  for  both  installing  of 
the  glass  using  epoxy  adhesive  and  for  low  temperature  testing.  This 
window  failed  during  cooldown.  The  frame  design  did  not  adequately 
limit  the  heat  transfer  from  the  outer  window  to  the  liquid  hydrogen 
heat  sink.  Consequently,  severe  frosting  of  the  window  occurred.  As 
a  result,  CVA  designed  a  new  window  configuration  which  was  reviewed 
by  NBS.  The  redesign  also  had  several  problem  areas,  in  particular, 
inadequate  fin  area  and  possible  weak  structure.  Detailed  comments 
on  this  new  design  were  submitted  in  a  letter  report  to  CVA. 

Assistance  to  several  other  CVA  departments  on  window  design 
for  low  temperature  applications  has  also  been  given  during  visits. 

The  discussions  with  CVA  personnel  on  this  general  subject  have  indi¬ 
cated  a  need  for  a  compilation  of  some  of  the  design  parameters 
necessary  for  a  successful  window.  An  informal  report,  NBS  Labora¬ 
tory  Note  61-2,  "Preliminary  Notes  on  Cryogenic  Windows"  has  been 
prepared  to  partially  fulfill  this  need.  A  window  to  demonstrate  the 
design  parameters  has  been  constructed  in  Boulder.  The  dimensions 
are  similar  to  those  required  for  the  Zero-G  test  vessel.  This  window 
will  be  tested  in  the  near  future.  An  additional  window,  using  the 
compressed  O-ring  seal  principle  reported  elsewhere  in  this  report, 
is  also  being  investigated.  A  test  apparatus,  to  determine  the  ability 
of  the  window  to  withstand  the  high  loading  forces  required,  is 
completed.  Results  of  the  tests  will  be  available  shortly. 
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Alt  hough  considerable  discussion  has  taken  place  with  the 
various  groups  at  CVA  concerned  with  window  designs,  a  number  of 
designs  have  recently  been  seen  in  the  CVA  shop  that  CEL.  feels  will 
fail. 

3.  2  Bearings  for  Zero-Gravity  Vent  Device 

The  possibility  of  using  ball  bearings  in  a  centrifugal  vent 
device  was  indicated  in  the  First  Quarterly  Report.  As  was  indicated 
in  that  report,  bearings  have  been  operated  successfully  in  a  cold 
gaseous  nitrogen  atmosphere  under  essentially  no  load.  It  is  necessary 
to  extend  the  testing  to  include  tests  under  conditions  similar  to  what 
is  expected  in  the  vent  device. 

To  evaluate  bearings  operating  in  a  hydrogen  atmosphere,  it 
was  necessary  to  relocate  the  NBS  test  apparatus  and  make  some 
modifications  to  the  existing  equipment.  The  tester  has  been  moved 
and  the  modifications  are  nearing  completion.  Preliminary  testing 
should  start  before  February  1. 

Actual  testing  will  be  concerned  with  finding  a  bearing  that  will 
operate  reliably  under  conditions  similar  to  what  are  expected  in  the 
venting  device.  The  tests  will  be  performed  with  the  bearings  running 
in  a  gaseous  hydrogen  atmosphere.  The  hydrogen  gas  will  be  cooled  in 
liquid  nitrogen  before  being  introduced  into  the  tester.  During  the 
tests  the  following  information  will  be  monitored:  speed,  load,  torque 
and  outer  ring  temperature  of  the  test  bearings. 

The  bearing  test  speed  will  be  approximately  10,  000  rpm,  the 
anticipated  operating  speed  of  the  vent  device.  Both  radial  and  thrust 
loads  will  be  encountered.  However,  only  thrust  loads  can  be  applied 
to  bearings  tested  in  the  existing  apparatus.  A  thrust  load  equivalent 
to  the  combined  radial  and  thrust  loads  that  are  expected  will  be  placed 
on  the  bearings.  Since  the  radial  load  is  quite  small  compared  to  the 
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thrust  load,  this  should  not  introduce  inconsistencies  between  test 
results  and  their  application.  Torque  will  be  measured  as  a  function 
of  time  with  a  variable  capacitor  torquemeter.  Temperature  will  also 
be  measured  as  a  function  of  time  and  coolant  gas  flow  rate  by  attach¬ 
ing  thermocouples  on  the  outer  rings  of  the  bearings. 

Possibilities  other  than  unlubricated  ball  bearings  appear 
feasible  for  the  centrifugal  vent  device.  Plans  are  presently  being 
made  to  extend  the  bearing  tests  to  include  tests  on  sleeve  type  journal 
bearings,  in  particular  bearings  made  of  carbon  compounds.  Prelim¬ 
inary  tests  will  be  conducted  submerged  in  liquid  nitrogen.  If  success¬ 
ful  bearings  are  established  under  these  conditions,  the  tests  will 
extend  to  operation  in  a  cold  gaseous  hydrogen  atmosphere. 

3.  3  Fuel  Line  Flange  Seal 

The  fuel  line  flange  currently  used  on  the  Centaur  employs  a 
pressure  actuated  "Teflon"  jacketed  metal  seal.  This  seal  has  failed 
to  perform  properly  in  this  application.  Recent  work*  performed  by 
NBS  reports  a  technique  which  has  made  excellent  seals  under  similar 
conditions. 

Accordingly,  an  investigation  to  study  the  use  of  this  seal  in 
the  Centaur  application  has  been  undertaken.  Two  test  flanges  (similar 
to  the  Centaur  flange)  were  constructed  as  illustrated  in  figure  3.  A 
test  container  was  also  constructed  to  allow  the  assembled  flange  to  be 
pressure  tested  at  liquid  nitrogen  temperature  (-320°F).  The  flange 
was  assembled  as  shown  in  figure  4  and  placed  inside  the  test  container. 
The  volume  around  the  assembly  was  evacuated  and  connected  to  a 
mass  spectrometer  helium  leak  detector.  The  volume  inside  the  flange 

♦Weitzel,  D.  H. ,  et  al,  "Elastomers  for  Static  Seals  at 
Cryogenic  Temperatures",  Paper  D-6,  Proceedings  of  the  I960  Cryo¬ 
genic  Engineering  Conf. ,  Boulder,  Colorado,  August  21-23,  I960. 
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Scale:  Full  Note:  S.  Fractional  tolerance*  1/64 

Material.*  321  SS.  2. Break  all  corners 

3.  Decimal  tolerance  *0.005 


Figure  3 


CENTAUR  FUEL  LINE  FLANGE 
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Figure  4 


FLANGE  ASSEMBLY 
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assembly  could  then  be  cooled  to  liquid  nitrogen  temperature  and 
pressurized  with  helium  gas.  Any  seal  leakage  would  immediately  be 
shown  by  the  leak  detector.  The  assembly  was  proof-tested  to  225 
psi  and  leak  checked  at  165  psi. 

Three  tests  have,  thus  far,  been  performed  with  this  apparatus. 
The  first  of  these  employed  a  neoprene  O-ring,  I.  D.  =  2.  614  in. , 
width  =  0.  070  in. ,  obtained  from  Plastic  and  Rubber  Products,  Inc. 
(Parco),  Los  Angeles,  California.  The  neoprene  was  designated  as 
compound  307-50.  The  seal  was  made  by  installing  the  O-ring,  as 
received  from  the  manufacturer,  between  the  two  clean  (degreased 
with  trichloroethylene)  flange  faces.  The  unit  was  assembled  by  evenly 
tightening  ten  1/4  inch  N.  F.  cap  screws  to  a  torque  of  120  in.  -lb. 

When  pressure  tested,  the  seal  was  found  to  be  leak  tight  on  the  rm.st 
sensitive  scale  of  the  mass  spectrometer  at  both  room  temperature 
and  liquid  nitrogen  temperature.  The  seal  assembly  was  temperature 
cycled  from  ambient  to  liquid  nitrogen  temperature  three  times,  and 
pressure  tested  after  each  cycle.  No  leaks  were  indicated  at  any 
time. 

Identical  tests  were  made  with  two  additional  O-rings: 

1.  Parco  Compound  921-65  ("Hypalon"),  I.  D.  =  2.  614  in. , 
width  =  0.  070  in.  Bolt  torque  =  70  in.  -lb. 

2.  Neoprene,  O.D.  =  2-3/4  in. ,  width  =  0.  139  in.  This 
standard  O-ring  was  a  rather  hard  compound  carried  in 
NBS-CEL  stock  for  general  high-vacuum  use.  Bolt 
torque  -  120  in.  -lb. 

Both  of  these  tests  were  as  pron  jing  as  the  first,  with  no  leaks 
indicated. 

It  is  interesting  to  note  that  the  flange  edges,  on  tests  one  and 
two,  were  touching  at  the  bolt  torque  indicated. 


-80- 


The  susceptibility  of  these  seals  to  vibrational  loads  was  not 
explored  in  these  preliminary  tests.  Previous  confined  seals  in 
similar  applications  have  demonstrated  ability  to  withstand  vibration. 

These  seals  appear  promising  for  the  Centaur  application. 

The  present  testing  was  done  at  liquid  nitrogen  temperatures  only. 

Past  experience  with  these  materials,  however,  indicates  that  service 
at  liquid  hydrogen  temperature  would  have  no  further  detrimental 
effects,  which  were  not  already  apparent  at  liquid  nitrogen  temperature. 

The  encouraging  results  warrant  further  investigation  of  the 
possible  application  to  the  fuel  line  flanges.  The  preliminary  results 
have  been  given  to  CVA  and  a  joint  program  established.  NBS  will 
perform  additional  thermal  and  pressure  tests  on  the  above  sealing 
method  and  on  other  more  conventional  commercial  seals.  CVA  will 
perform  vibration  tests  on  the  successful  seals.  Physical  retainment 
of  the  seals  for  field  installation  will  be  considered  by  NBS  and  CVA. 

3.  4  Ground  Support  Equipment 

Several  general  discussions  were  held  with  CVA  personnel 
concerning  ground  support  equipment.  These  had  particular  applica¬ 
tion  to  facilities  at  Point  Loma  and  Sycamore  Test  Sites.  The  subjects 
of  transfer  lines,  storage  tanks,  high  vacuum  insulation  and  powder- 
vacuum  insulation  were  considered. 

A  contamination  problem  has  been  encountered  involving  the 
pressurizing  gas  used  for  fuel  transfer.  The  helium  gas  used  for  this 
purpose  contains  several  hundred  parts  per  million  of  oxygen  and 
nitrogen  in  the  proportions  of  air.  The  potential  hazards  of  the 
accumulated  air  and  methods  of  avoiding  the  initial  contamination  were 
discussed.  Based  on  past  experience  and  experimental  programs 
concerned  with  the  potential  hazards  of  solid  air  in  liquid  hydrogen, 
it  was  concluded  that  although  no  particular  hazard  existed,  steps 
should  be  taken  to  prevent  the  initial  contamination  of  the  helium. 
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3.  5  Propellant  Tank  Insulation 

Insulation  for  the  propellant  tank  has  application  to  the  Propel¬ 
lant  Test  Vehicle  (PTV)  and  the  Centaur  Vehicle.  The  PTV  insulation 
is  quite  heavy  compared  to  the  Centaur  Vehicle  in  order  to  simulate 
the  heat  transfer  encountered  in  space.  It  is  of  immediate  concern 
•  to  obtain  a  successful  insulation  for  such  application.  The  Centaur 

insulation  is  divided  into  that  which  is  jetisoned  during  flight  and  that 
«  which  remains  aboard  on  the  fuel  tank  forward  bulkhead.  Success  of 

the  insulation  is  of  utmost  importance. 

During  the  reporting  period,  sprayed  on  foam  was  applied  to 
two  PTV  tanks.  The  PTV  located  at  Sycamore  has  been  cooled  on 
several  occasions.  On  the  first  cooldown  the  insulation  exhibited 
typical  failures  --  circumferential  cracks,  radial  cracks  and  poor 
adhesion.  On  subsequent  cooldowns,  the  initial  failures  propagated. 
The  effectiveness  of  the  insulation  is  quite  poor. 

Although  no  immediate  solution  to  the  problem  is  apparent,  it 
is  felt  that  further  study  should  not  be  de -emphasized.  Suggestions  on 
a  possible  program  of  research  and  development  have  been  prepared. 
Further  study  should  include  not  only  attempts  to  obtain  an  immediate 
solution,  but  also  a  long-range  development  program. 

4.  A  Compilation  of  Thermophysical  Properties  of  Cryogenic 
Materials 

R.  B.  Stewart  amd  V.  J.  Johnson 

j 

The  task  of  compiling  data  of  cryogenic  materials  was  started 
at  CEL.  in  January  1958.  This  early  effort  culminated  with  the  publi¬ 
cation  of  the  preliminary  report,  "A  Compendium  of  the  Properties  of 
Materials  at  Low  Temperature",  Phase  I,  in  December  1959.  The 
scope  of  this  activity  was  later  increased  in  January  1959  to  include 
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seven  additional  properties  which  was  the  beginning  of  the  work  now 
being  carried  on  in  this  project. 

Extensive  bibliographies  have  been  prepared  for  this  second 
phase  of  this  activity  on  the  following  subjects: 

1.  Compressibility  Factor  of  Cryogenic  Fluids. 

2.  Compressibility  and  Compressibility  Coefficients  of 
Cryogenic  Fluids. 

3.  Entropy  Data  for  Cryogenic  Fluids. 

4.  Velocity  of  Sound  for  Cryogenic  Fluids. 

5.  Solubility  of  2- Component  Mixtures  of  Cryogenic  Fluids. 

6.  Electrical  Resistivity  (and  Conductivity)  of  Materials  at 
Low  Temperatures. 

7.  Ferromagnetic  Properties  of  Materials  at  Low  Temperatures. 

(The  cryogenic  fluids  of  primary  interest  are:  Helium,  Hydro¬ 
gen,  Neon,  Nitrogen,  Oxygen,  Air,  Carbon  Monoxide,  Fluorine, 

Argon  and  Methane. ) 

A  task  of  assembling  thermal  conductivity  integrals  of  cryogenic 
materials  was  completed  from  work  done  in  Phase  I.  (The  thermal 
conductivity  integral  is  a  cumulative  value  of  thermal  conductivity 
from  the  datum  temperature  rather  than  an  instantaneous  value. )  This 
material  is  compiled  in  tables  for  44  pure  metallic  substances,  36  non- 
ferrous  alloys,  9  ferrous  alloys  and  4  glasses  and  plastics. 

The  compilation  and  correlation  of  data  and  the  presentation  of 
the  results  on  completely  documented  data  sheets  in  the  current 
established  program  is  progressing.  These  data  sheets  usually  contain 
graphical  presentation  of  the  data  in  a  form  that  is  particularly  useful 
for  engineering  information.  In  addition,  the  data  sheet  includes  a 
documentation  of  the  data  sources,  information  on  methods  of  analysis 
of  the  data  and  basis  for  selection  of  the  data  used  together  with 
comparisons  of  alternate  sources  when  available.  Tables  of  selected 
values  are  also  included. 
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Data  sheets  are  currently  being  prepared  for  the  subjects 
included  in  the  above  bibliographies.  These  data  sheets  are  in  various 
stages  of  completion  for  the  materials  currently  being  considered. 
Many  of  the  data  sheets,  in  particular  on  compressibility  factor, 
velocity  of  sound,  solubility  and  electrical  resistivity  are  in  the  final 
stages  of  review  for  issuance. 

The  following  is  an  indication  of  the  current  status  of  tasks 
worked  on  during  the  reporting  period,  together  with  an  estimate  of 
their  completeness: 

1.  Compressibility  Factor  for  Cryogenic  Fluids  (Z  =  PV/RT) 
(data  sheets  for  hydrogen,  nitrogen,  methane  and  air  were 
completed  previously) 

Argon  (1  to  5000  atm.,  120  to  300°K)(70%  completed) 

Helium  (1  to  100  atm.,  20  to  300°K)(50%  completed) 

Neon  (20  to  90  atm.,  55  to  300°K)(100%  completed) 

2.  Entropy  for  Cryogenic  Fluids  (T-S  diagram  with  constant 
property  lines  for  pressure,  volume,  enthalpy) 

Neon  (20  to  90  atm. ,  55  to  300®K)(60%  completed) 

Helium  (1  to  100  atm.,  20  to  300°K)(30%  completed) 
Nitrogen  (1  to  3000  atm.,  90  to  300°K)(20%  completed) 

3.  Solubility  and  Phase  Equilibria  of  2-Component  Mixtures 
(solid,  liquid,  gas  phases).  About  28  data  sheets  are  in 
various  states  of  preparation.  (This  task  is  now  about 
80%  completed. ) 

4.  Electrical  Resistivity  of  the  Pure  Metals.  Data  sheets  for 
21  metals  have  been  completed  for  review.  Data  sheets 
for  25  metals  are  in  progress  and  are  about  80%  completed. 

In  addition  to  continuing  the  tasks  now  in  progress  to  the 
completion  of  several  data  sheets  in  the  next  reporting  period,  it  is 
anticipated  that  active  preparation  of  data  sheets  will  be  undertaken 
for  the  subjects  for  which  bibliographies  are  now  available.  Additional 
materials  and  an  increased  number  of  fluids  will  also  be  considered, 
providing  adequate  staffing  can  be  made  available. 
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*  ELECTRICITY.  Resistance  and  Reactance.  Electroelm mistry.  Electrical  Instruments.  Magnetic  Measurements. 

Dielectrics. 

METROLOGY.  Photometry  and  Colorimetry.  Refraetometry.  Photographic  Research.  Length.  Engineering 
^  Metrology.  Mass  and  Scale.  Volumetry  and  Densimetry. 

HEAT.  Temperature  Physics.  Heat  Measurements.  Cryogenic  Physics.  Rheology.  Molecular  Kinetics.  Free 
Radicals  Research.  Equation  of  State.  Statistical  Physics.  Molecular  Spectroscopy. 

RADIATION  PHYSICS.  X-Ray.  Radioactivity.  Radiation  Theory.  High  Energy  Radiation.  Radiological 
Equipment.  Nucleonic  Instrumentation.  Neutron  Physics. 

CHEMISTRY.  Surface  Chemistry.  Organic  Chemistry.  Analytical  Chemistry.  Inorganic  Chemistry.  Electro¬ 
deposition.  Molecular  Structure  and  Properties  of  Gases.  Physical  Chemistry.  Thermochemistry.  Spectrochem¬ 
istry.  Pure  Substances. 

MECHANICS.  Sound.  Pressure  and  Vacuum.  Fluid  Mechanics.  Engineering  Mechanics.  Combustion  Controls. 
ORGANIC  AND  FIBROUS  MATERIALS.  Rubber.  Textiles.  Paper.  Leather.  Testing  and  Specifications. 
Polymer  Structure.  Plastics.  Dental  Research. 

METALLURGY.  Thermal  Metallurgy.  Chemical  Metallurgy.  Mechanical  Metallurgy.  Corrosion.  Metal  Physics. 
MINERAL  PRODUCTS.  Engineering  Ceramics.  Glass.  Refractories.  Enameled  Metals.  Constitution  and 
Microstructure. 

BUILDING  RESEARCH.  Structural  Engineering.  Fire  Research.  Mechanical  Systems.  Organic  Building  Ma¬ 
terials.  Codes  and  Safety  Standards.  Heat  Transfer.  Inorganic  Building  Materials. 

APPLIED  MATHEMATICS.  Numerical  Analysis.  Computation.  Statistical  Engineering.  Mathematical 
Physics. 

DATA  PROCESSING  SYSTEMS.  Components  and  Techniques.  Digital  Circuitry.  Digital  Systems.  Analog 
Systems.  Applications  Engineering. 

ATOMIC  PHYSICS.  Spectroscopy.  Radiometry.  Mass  Spectrometry.  Solid  State  Physics.  Electron  Physics. 
Atomic  Physics. 

INSTRUMENTATION.  Engineering  Electronics.  Electron  Devices.  Electronic  Instrumentation.  Mechanical 
Instruments.  Basic  Instrumentation. 
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CRYOGENIC  ENGINEERING.  Cryogenic  Equipment.  Ciyogenic  Processes.  Properties  of  Materials.  Gas 
j  Liquefaction. 

IONOSPHERE  RESEARCH  AND  PROPAGATION.  Low  Frequency  and  Very  Low  Frequency  Research.  Ion¬ 
osphere  Research.  Prediction  Services.  Sun-Earth  Relationships.  Field  Engineering.  Radio  Warning  Services. 
RADIO  PROPAGATION  ENGINEERING.  Data  Reduction  Instrumentation.  Radio  Noise.  Tropospheric  Meas- 
r  urements.  Tropospheric  Analysis.  Propagation-Terrain  Effects.  Radio-Meteorology.  Lower  Atmosphere  Physics. 

RADIO  STANDARDS.  High  frequency  Electrical  Standards.  Radio  Broadcast  Service.  Radio  and  Microwave 
Materials.  Atomic  Frequency  and  Time  Standards.  Electronic  Calibration  Center.  Millimeter-Wave  Research. 
Microwave  Circuit  Standards. 

RADIO  SYSTEMS.  High  Frequency  and  Very  High  Frequency  Research.  Modulation  Research.  Antenna  Re¬ 
search.  Navigation  Systems.  Space  Telecommunications. 

UPPER  ATMOSPHERE  AND  SPACE  PHYSICS.  Upper  Atmosphere  and  Plasma  Physics.  Ionosphere  and 
Exosphere  Scatter.  Airglow  and  Aurora.  Ionospheric  Radio  Astronomy. 


